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ABSTRACT 

In this paper we propose that the compressive tidal field in the centers of flat-core 
early type galaxies and ultraluminous galaxies compresses molecular clouds producing 
dense gas obseved in the centers of these galaxies. The effect of galactic tidal fields 
is usually considered disruptive in the literature. However, for some galaxies, the 
mass profile fiattens towards the center and the resulting galactic tidal field is not 
disruptive but instead it is compressive within the fiat-core region. We have used the 
virial theorem to determine the minimum density of a molecular cloud to be stable 
and gravitationally bound within the tidally compressive region of a galaxy. We have 
applied the mechanism to determine the mean molecular cloud densities in the centers 
of a sample of fiat-core, early-type galaxies and ultraluminous galaxies. 

For early-type galaxies with a core-type luminosity profile, the tidal field of the 
galaxy is compressive within half the core radius. We have calculated the mean gas 
densities for molecular gas in a sample of early type galaxies which have already been 
detected in CO emission and we obtain mean densities of (n) ~ 10^ — lO^cm^^ within 
the central lOOpc radius. We also use our model to calculate the molecular cloud 
densities in the inner few hundred parsecs of a sample of ultraluminous galaxies. From 
the observed rotation curves of these galaxies we show that they have a compressive 
core within their nuclear region. Our model predicts minimum molecular gas densities 
in the range 10^ — lO^cm"^ in the nuclear gas disks; the smaller values are applicable 
typically for galaxies with larger core radii. The resulting density values agree well with 
the observed range. Also, for large core radii, even fairly low density gas (~ lO^cm"^) 
can remain bound and stable close to the galactic center. 



Subject headings: galaxies : elliptical and lenticular - galaxies : kinematics and 
dynamics - galaxies: nuclei - ISM : clouds 
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1. INTRODUCTION 

Over the last two decades observations have revealed the presence of significant amounts of 
cold, neutral gas in early type galaxies (Henkel &; Wiklind 1997; Knapp 1998). The term "early 
type galaxies" includes lenticular or SO galaxies, dwarf elliptical (dE) and elliptical (E) galaxies. 
Many of these galaxies have prominent dust lanes and extended distributions of neutral hydrogen 
gas (HI). A significant fraction of these galaxies also contain molecular hydrogen gas (H2). Unlike 
the HI gas which is usually distributed in the outer parts of the galaxy, the molecular gas is nearly 
always centrally concentrated. The H2 gas in these galaxies has been detected and traced using 
the CO(l-O) and C0(2-l) lines, indicating that the gas may be fairly dense lO^cm"^). In a 
few cases high density molecular gas n(H2) > lO^cm"^ has been observed using the high density 
tracers HCN and HCO"*". The amount of molecular gas observed in early type galaxies varies 
between IO^Mq to IO^Mq and is less than in spirals of the same blue luminosity (Knapp 1998). 
While the amount of cold gas is small compared to that in spiral galaxies, it is still not insignificant 
and in some cases also forms the site of ongoing star formation. The origin and spatial distribution 
of this dense gas has not been explained so far. 

Ultraluminous galaxies are very bright starburst galaxies having far-IR luminosities of the 
order L > IO^^l^ (Sanders et al. 1988). They are formed from the merger of two gas rich spiral 
galaxies. The merger is expected to result in large amounts of disk molecular gas falling into the 
center of the newly formed merger remnant, which subsequently results in enhanced star- formation 
or super-starburst at the center. The total central gas mass is ~ 10^ — lO^'^M©, and is distributed 
over the central few hundred pc. 

It may seem strange at first that wc are studying seemingly different kinds of galaxies (early 
type galaxies and ultraluminous galaxies) together in this paper. However, as we argue next, these 
are close in terms of their dynamics. Mergers of two spirals are known to result in an ultraluminous 
galaxy, as seen from observations (Sanders et al. 1988), and as explained theoretically (e.g.. 
Jog & Das 1992). A merging pair is believed to be a precursor of an elliptical galaxy, as shown 
observationally (e.g., Schweizer 1999), and in fact the K-band brightness profiles of about 50 % of 
the ultraluminous galaxies display a r^/^ distribution typical of an elliptical galaxy (Sanders & 
Mirabel 1996). Numerical simulations have shown that a merger of two spirals can result result 
in an elliptical galaxy (e.g., Barnes & Hernquist 1992). This scenario is supported by the fact 
that the observed gas surface density in the centers of ultraluminous galaxies is 5OOM0pc~^ and 
is comparable with the central stellar density in ellipticals (Kormendy k, Sanders 1992; Solomon 
1997). Therefore, wc expect that the mass distribution and hence the dynamics of the centres of 
these galaxies would be similar. Also, the gas number densitites would be affected in a similar way 
by the galactic tidal field. Hence we study the origin of the dense gas in the centres of early type 
and ultraluminous galaxies together in this paper. We do not, however, study the relevance of the 
compressive tidal field, which we discuss next, for the evolution from a merging pair of spirals into 
an elliptical galaxy. 
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The galactic tidal field is usually thought to be disruptive in the centers of galaxies so that 
only dense gas would survive the strong shearing due to a central tidal field. However, if the 
potential flattens instead of peaking towards the center of the galaxy, the tidal field can reverse 
sign and become compressive. The tidal field is compressive when the background gravitational 
field of the galactic disk enhances the self-gravitational field of the cloud (Binney Sz Tremaine 
1987, Chap. 7). We show this situation is physically possible in core regions of some galaxies - 
such as the flat-core early-type and ultraluminous galaxies. This is a new application of the idea 
of compressive tidal fields, namely, in the astrophysical context of galaxy centers. In the past, the 
effect of the compressive nature of tidal fields has been studied in the context of disk shocking 
of globular clusters (Ostriker, Spitzer, & Chevalier 1972), and the tidal effect on a galaxy in the 
core of a cluster of galaxies (Valluri 1993). The other new feature of our study is the application 
of compressive fields to a compressible fluid, namely gas, and our calculation of the steady-state 
density of a bound, virialized cloud in the presence of this field. 

In §2 and §3, we determine the minimum density of a stable cloud in a compressive tidal 
field. We have applied these results (§4) to the centers of a few fiat-core, early-type galaxies in 
which molecular gas has been detected and also to a sample of ultraluminous galaxies. We have 
compared our results with the observed gas densities. We thus show that it is the compressive 
nature of the galactic tidal fields in the centres of these galaxies that naturally leads to the origin 
of the dense molecular gas observed in these galaxies. 



2. EQUILIBRIUM OF A CLOUD IN THE TIDAL FIELD OF A GALAXY 

Molecular clouds are generally bound by self-gravity, hence the compressive galactic tidal field 
plays an important role in the stability of a cloud. In this section we use the galactic tidal field in 
a spherical potential to derive the net stability criterion for a cloud. 



2.1. Tidal Field 



The components of the tidal field in a spherical potential has been derived in some detail by 
Cowsik &: Ghosh (1987), and Valluri (1993) who studied the effect of the compressive tidal field of 
a cluster potential on a galaxy. We summarize their results briefiy here before applying them to a 



molecular cloud. For a potential $(Xi) the tidal field components at xi, where Xi 
given by. 



Xi — Xi, 



dXidXi 



are 



(1) 



For a spherical potential $(R) for which the cloud lies along the X axis so that X=R, Y=Z=0, the 
components of the tidal field at a position (x,y,z) in the coordinate system centered on the cloud. 
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are given by, 



X 



'dR 



(2) 

(3) 
(4) 



Let {p) be the mean density of tlie mass witliin radiTis R in the galaxy so that [p) = ^^f^, where 
M(R) is the mass of the galaxy within radius R. Then the tidal field components become, 



-7rG(2(p) - ?,p)x 



TrG{p)y 



(5) 
(6) 
(7) 



T, = --7rG{p)z 

Thus, the tidal field is always compressive along the directions perpendicular to the radius 
vector for a spherical potential. It can be either compressive or disruptive along the radius 
vector depending on the density distribution in the galaxy. For p{R) > the tidal field is 
compressive. The nature of the tidal field thus depends on the density profile of the system. For a 
centrally peaked density profile, the tidal field will always be disruptive along the radial direction. 
However, if the density profile of the galaxy flattens towards the center, then the tidal field can 
be compressive within the flat core region of the galaxy. The tidal field is compressive in the 
centers of clusters of galaxies and is important in the evolution of galaxies in the cluster (Valluri 
1993). We show that it is also compressive in the centers of flat core early type galaxies and some 
ultraluminous galaxies; this is derived using equation (2) and is discussed in Sections 3 and 4 of 
this paper. Equations (5)-(7) are shown to illustrate the physics of compressive tidal fields. For 
these galaxies the tidal field changes direction from being a disruptive force to a compressive force 
in the center of the galaxy. We show that this is an important factor in deciding the steady state 
density of a cloud in the galaxy center. 

When the tidal field is compressive along the line of action of the force, the cloud will feel an 
overall compressive force. In such cases the tidal field can never disrupt the cloud, instead the 
cloud will be compressed by the galactic field. For a cloud in equilibrium in a compressive tidal 
field, the internal pressure maintains equilibrium against both the self gravity of the cloud and the 
compressive tidal field of the galaxy. 
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2.2. Virial Equilibrium of a Cloud 

For a system in virial equilibrium, we have the relation 

-2W=EFi.ri (8) 

i=l 

where (K) is the time averaged kinetic energy of the system, Fi is the force on the i*^ particle 
and r; is the position coordinate of the i*^ particle wrt the cloud center. For a molecular cloud in 
virial equilibrium, the forces involved will be the internal self gravity of the cloud and the galactic 
tidal field acting on the cloud. As mentioned previously, the tidal field is anisotropic. However for 
simplicity we will consider only the radial component and not the transverse components of the 
tidal field in the deriving the equilibrium density of a cloud. Strictly speaking we should include 
all three components of the tidal field which would lead to an ellipsoidal cloud. However, in a real 
cloud, the compression would tend to become isotropic soon (see Section 5). Hence, for simplicity, 
we assume that the cloud is spherical and we include only the radial tidal field for estimating the 
equilibrium cloud density. If we convert the right hand side of the equation (8) into an integration, 
we obtain the following relation, 

2 rGMJr)dMJr) _ /" o,,. 
- Mcv^ = - ^ +To r^dMc (9) 

where Mc is the mass of the cloud, v is the 3-D velocity dispersion in the cloud and Tq is the tidal 
field per unit distance across the cloud (i.e. Tq = -0). We assume that the density within the 
cloud is inversely proportional to radius as observed (Sanders, Scoville, &; Solomon 1985), so that 
Pc(r) = where po is the density of the cloud at the cloud edge, i.e. at r = Rp. Then the 

cloud has a mass profile given by Mc(r) = 27rRcPor^- Thus, the mean density of the cloud, {pc), is 
obtained to be (pc) = |po- The equation of cloud equilibrium thus reduces to. 

This has been derived in the galactic rest frame. In the Appendix A, we give the derivation 
in the rest frame of the cloud for the sake of completeness, and argue that for considering cloud 
stability it is appropriate to use the galactic frame equation as derived here. 



2.3. Minimum Density for Stability 

For a tidally disruptive field, the tidal field coefficient Tq is positive. Since v^ > 0, we obtain 
the familiar relation for tidal stability of a molecular cloud, namely 

2GMc 1 9 
3^ > 2^°^ 
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Writing this condition in terms of (pc), the mean density of the cloud, we get 

te) > -i^n (12) 

Thus, the minimum density for a cloud to be stable against the disruptive tidal field is 

= -i^n (13) 

For a tidally compressive field Tq = — | Tq |< 0, hence v^ > does not give us any new information 
except that the cloud can never be disrupted by the tidal field. So to obtain the density of a cloud 
in the compressive tidal field we impose the condition that the cloud be gravitationally bound to 
begin with, i.e. 



Using equation (10) we get 



= 16^ I ^» I 

which is formally the same relation as that for a disruptive tidal field the only difference being that 
we take the magnitude of the tidal force coefficient. Physically, however they have very different 
significance. For a disruptive field, the minimun density (pm) is the density below which a cloud 
will be disrupted by the galactic tidal field. For a compressive tidal field it denotes the minimum 
density for a gravitationally bound cloud to exist within that region. Since the cloud is assumed 
to be gravitationally bound to begin with, an additional inward/compressive field increases the 
velocity dispersion required to satisfy the virial equilibrium (eq.[10]), and hence the critical density 
required for stability in presence of a compressive force is actullay higher for a virialized cloud. 
For lower densities, the gas within that region will be diffuse and gravitationally unbound. Such 
gas will not lead to star formation. 

Note that the critical number density of hydrogen molecules (in cm~^) corresponding to 

the mass density (pm) (eq.[15]), assuming a 10 % gas number fraction in helium, is given by 
(Pm)/(4.68 X 10^^^). We will use this conversion in the next two sections. To calculate and plot 
tidal fields, we use mass, length and time units to be M©, pc and lO^year respectively. This means 
that the units of Tq are ( iqt^) ■ 



3. APPLICATION TO MOLECULAR GAS IN EARLY TYPE GALAXIES 

In this section we apply the above criterion to predict the central molecular gas densities in 
six early type galaxies. As mentioned later in this section, we need parameters that can be derived 
from the luminosity profile. Unfortunately, only a few of the galaxies whose surface brightness 
has been studied have also been detected in CO. Hence we could apply the results of the previous 
section to only a small sample of early type galaxies. 
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3.1. Spherical Stellar Potential 

Trcmainc ct al. (1994) have derived a one parameter set of models for spherical stellar 
systems, these give self-consistent potential-density pairs. Their model can be applied to the center 
of early type galaxies and spiral bulges. Prom their model they derive the density distribution, 
surface brightness profile and line of sight velocity dispersion. For all the models, at large distances 
the density distribution tends to /j(r) oc r~^ and the surface brightness tends to the familiar r^/^ 
form. The family of potentials, called rj models are given by. 



$(r ) — 



(1 + r')''-^ 



(16) 



where r' = R/Rq, < < 3, Mg is the mass of the galaxy and Rq is the core radius. Many of the 

properties of the i] models have also been derived by Dehnen (1993). The rj = 1 model corresponds 
to the Jaffe (1983) model while the i] = 2 model corresponds to the Hernquist (1990) model. Prom 
equation (16), Tq, the tidal field per unit distance is found to be, 

^0 = -^^(I^(.-2-2/) (17) 

For < Tj < 2 the density diverges towards the center of the galaxy and the tidal field is always 
disruptive. In contrast, for 2 < < 3 the model has a more flat density core and the tidal field 
is compressive within the core. However, even here, the field is compressive only for distances 
R < ^Rq. 

Thus to determine the critical density in the centers of these galaxies we need to determine 
two parameters; one is Mg, the mass of the galaxy, which can be determined from the mass to 
luminosity ratio (M/Lv) and the other is the observed core or break radius Rq. The break radius 
can be determined from the surface brightness profile of the galaxy. For a flat core profile the 
break radius is very sharp but for a power law profile, which is rising towards the center, Rq is not 
so sharply defined. 



3.2. Results for Minimum Gas Densities 

There are now a very large number of early type galaxies with detected H2 gas (Knapp 1998; 
Henkcl & Wiklind, 1997). In fact, after looking at the literature, we found that there are over 
fifty elliptical and SO galaxies with detected H2 gas. Many show both CO(l-O) and C0(2-l) lines 
towards the center indicating that there may be fairly dense gas in these galaxies. The fact that 
dense gas does exist in many early type galaxies indicates that there must be some stabilising 
mechanism resulting in fairly long lived dense gas in the centers of these galaxies. 

In this section we use the results of the previous section to explain the observed molecular 
gas densities in a sample of elliptical and SO galaxies. To determine (pm) we need both the 
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galaxy mass and the break radius to apply the analytical models of Tremaine et al. (1994). 

These parameters have been determined for a large sample of early type galaxies by Faber et al. 
(1997). They used very high resolution HST observations to determine the central parameters of 
61 early type galaxies. We compared their sample with CO observations of elliptical galaxies in 
the literature. We found that molecular gas has been detected in the center of only a few of the 
galaxies in the Faber et al. (1997) sample. We have used the data in column 3 and 5 of Table 3 in 
Faber et al. (1997) to determine the mass of the galaxy Mg and we used the data in Table 2 of 
their paper to determine the break or core radius Rq of the galaxy. The results are shown in Table 
1 where we have listed the galaxies, the parameters Rq and Mg and the density of molecular gas 
predicted from equation (15) of the previous section. Figure 1 illustrates how the tidal field varies 
in a power law galaxy (NGC 4697) and in a core profile galaxy (NGC 4649). For the latter, the 
tidal field becomes negative i.e. compressive within a certain radius in the galaxy, while for NGC 
4649 the tidal field continues to rise towards the center. Figure 2 shows the variation of miminum 
density in the disruptive and compressive tidal regions of NGC 4472. The dashed line demarcates 
the regions of the compressive and disruptive tidal field. In the following paragraphs, we have 
compared our results with the CO observations of the sample of galaxies taken from Faber et al. 
(1997). 

NGC 4594 and NGC 4697 : Both these galaxies have power law profiles and thus the tidal field 
is always disruptive in both of them. NGC 4594 (also called 'Sombrero Galaxy') has a very steep 
brightness profile in the center and also has an active nucleus (Bajaja et al., 1988). Bajaja et 
al. (1991) observed this galaxy in the CO (1-0) and (2-1) line and though they detected CO in 
the disk of the galaxy, they could not detect any CO in the center. The reason for this could be 
that the tidal field is very high in the inner few hundred pc and the gas cannot survive there. 
NGC 4697 on the other hand has a moderately steep luminosity profile with a break radius of 
~ 240pc (Faber et al. 1997). Using 77 = 1.5 in equation (17) we get the tidally limited density to 
be (n) >~ 5 x lO^cm"^ at 1 kpc and it rises to (n) ~ lO^cm"^ at 200pc. Sofue and Wakamatsu 
(1993) have detected CO emission from the central 1 kpc of this galaxy. Since the mass profile 
does not rise very steeply in the center of this galaxy, the tidal field is not very strong and so gas 
can exist in the center of NGC 4697 even though it has a power law profile like NGC 4594. We do 
not address the origin of high density gas observed in some power-law galaxies. 

NGC 4472 : This galaxy has a flat core luminosity profllc and hence the tidal field is compressive 
in the center. The break radius is ~180pc which means that gas will feel the compressive tidal 
field within the inner lOOpc where densities of at least (n) ~ lO^cm"^ are predicted by equation 
(15). CO emission detection indicates that there is ^ IO^Mq of molecular gas in this galaxy 
(Huchtmeier et al. 1994; Lees et al. 1991). However, Huchtmeier et al. (1988) could only detect 
CO from an off center position and not at the center. One reason for this could be that their beam 
was very large (~ 2kpc). Higher resolution observations might lead to detection of high density 
gas in the center of this galaxy. 

NGC 4649 : This galaxy is a giant elliptical belonging to the Virgo cluster. It has a fairly large 
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core which is compressive at radii less than 130pc. The model predicts gas densities of at least 
< n >~ 5 X lO^cm^'^ in the inner 130pc. C0(2-l) and CO(l-O) emission has been detected from 
this galaxy (Lees et al. 1991; Sage & Wrobel, 1989). 

NGC 1400 & NGC 1316 : These galaxies both have small cores and the model predicts stable gas 
of densities of at least (n) >'-^ lO^cm"^ to lO^cm^^ in the central 10 to 20pc of these galaxies. If 
the galaxy had a power law profile the gas would probably not be stable so close to the center. 
CO emission has been detected from NGC 1400 by Lees et al. (1991) and in NGC 1316 by Sage 
and Galletta (1993). 

NGC 1600, NGC 3379 & NGC 4486 : These galaxies have not been detected in CO, only upper 
limits to their molecular gas masses have been obtained (Sofue k. Wakamatsu, 1993). All three 
galaxies are ellipticals with compressive cores and hence dense gas can exist in the inner lOOpc of 
their cores. 



4. APPLICATION TO MOLECULAR GAS IN ULTRALUMINOUS GALAXIES 

Ultraluminous galaxies are merger remnants with large amounts of molecular gas concentrated 
in their centers. Over the past few years observations have revealed that the nuclear molecular gas 
is located in a central disk of size smaller than a few hundred pc (Scoville et al. 1991; Downes, 
Solomon k, Radford, 1993). Downes & Solomon (1998) observed with high resolution a sample 
of ultraluminous galaxies and showed that most of the CO flux from the center comes from 
subthermally excited CO emission. They concluded that most of the gas is of density ~ lO^cm"^ 
and only --^25% of it is high density gas (> lO^cm"^). We use our model to determine the 
minimum density of bound, virialized molecular gas clouds for their sample of galaxies using their 
parameters for the rotation curves, and show that our resulting values compare well with their 
observations. Despite the high observed turbulent velocities (see Table 5, Downes & Solomon, 
1998), the assumption of an individual cloud being bound may still be applicable if the beam 
covers a number of clouds and hence a large fraction of the linewidth is due to cloud-cloud velocity 
dispersion. Also, the volume filling factor of the central molecular gas is ~ 1. Nevertheless, we 
have treated the central gas as being discrete bound clouds for simplicity. 



Downes &; Solomon (1998) have observed the rotation curves of the molecular gas disks in 
a sample of ten ultraluminous galaxies. They found that the rotation curves rose to a steady 
velocity Vq at a radius and then became flat at the steady value. They modelled the rotation 
curves with the simple expression. 



4.1. 



Tidal Field from the Nuclear Rotation Curve 




(18) 
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where /3 = 1 for R < Rd and /3 = for R > Rd upto the maximum radius at which CO was 
observed. 

For the power law rotation curve (for R < R^) , the corresponding potential can be calculated 
to be: 

Hence, Tq, the tidal field along the radial direction is obtained to be: 

To = -(2/3-1)^ (20) 

Note that Tq < 0, that is the radial component of the tidal field is compressive, when /? > 1/2. 
Therefore, for /? = 1 as seen in the inner regions of ultraluminous galaxies, the galactic tidal field 
is compressive. For the region outside of Rj where the rotation curve is flat, the galactic tidal field 
is not compressive, instead it is disruptive. 

In Figure 3, we show the plot of the tidal field per unit mass versus the galactic radial distance 
for Arp 193 for the data from Downes &; Solomon (1998). The tidal field is compressive inside of 
220 pc. Note that sharp transition from disruptive to compressive region is an artifact due to the 
rotation curve model obtained by Downes & Solomon (1998), the actual transition is probably 
smoother. 



4.2. Results for Minimum Gas Densities 

We used the equation (15) and equation (20) to determine the minimum density for a bound, 
virialized cloud to exist in the compressive tidal field of the ultraluminous galaxies observed by 
Downes & Solomon (1998). Our results are shown in Table 2 where we have listed the galaxies, 
their rotation curve turnover velocity, turnover radius and the minimum gas density < Pm >• The 
particle number density (n) lies in the range of 10^ to 10"^ cm . This density is the minimum 
density of a gravitationally bound cloud to remain stable in the center of the galaxy. This range 
agrees very well with the mean density of the main mass component of the molecular gas observed 
in these galaxies (see Table 6, Downes & Solomon 1998). Since we have assumed discrete clouds 
with a volume filling factor < 1 (§4), our values are slightly larger than the observed values, as 
expected, since the latter have been derived assuming a uniformly filled disk. 

The peak density of the molecular gas is observed to be near the turnover point in the 
rotation curve whereas our model predicts a constant high density inside this point. Also there 
is substantial gas observed beyond the turnover point, but both these could arise because of the 
sharp turnover in the model rotation curve assumed by Downes &; Solomon (1998). 

Note that the very dense gas with n > lO^cm"^ which constitutes about 25% of the gas mass 
in the extreme starburst regions such as Arp 220 West and East could also arise due to tidal 
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compression (see Table 12, Downes k, Solomon 1998 and our Table 2). This high density gas is 
the site of HCN emission, and contains a large fraction of the star formation. Since the internal 
velocity dispersion increases with density for a bound virialized cloud, the observed velocity 
dispersion in a given region sets an upper limit on the high density that can result from this 
mechanism. The observed dispersion in a region could be due to the superposition of a number 
of moderate density clouds ((n) ~ 10^ — lO^cm"^) in the beam, as we argued in §4, or it could be 
due to a few higher density clouds. We do not discuss the issues of the fraction of the gas in the 
dense component, and the origin of the high central turbulent velocities, in this paper. 

It is interesting that even low to moderate density gas (n) ~ 10^ — lO^cm"^ in some cases can 
survive in the center and will not be disrupted by the galactic tidal shear. Downes & Solomon 
(1998) suggest that most of the moderate density gas (10^cm~^) that they observe in these galaxies 
is diffuse and unbound as it may not be able to withstand the tidal shear of the nuclear region 
of the galaxy. However as we have shown here, gas can survive near the center as bound clouds 
because the tidal field is compressive and not disruptive. Hence the gas that they have termed as 
the 'diffuse medium' could well be bound clouds. 

Since these clouds are bound, they are more susceptible to the onset of star formation. Thus, 
our conclusion that these clouds are bound helps the scenario proposed by Downes &; Solomon 
(1998) for starbursts powering the ultraluminous galaxies. The starbursts could be triggered due 
to large-scale gravitational star-gas instability (Jog 1996), as they suggest, and/or it could arise 
due to the collapse of the individual, bound clouds, or due to collisions between the bound clouds. 
Hence the formation of dense gas ((n) > 10^ — lO^cm"^) in a compressive tidal field is important 
for the starburst and hence the evolution of ultraluminous galaxies. 

5. DISCUSSION 

1. Normal spirals also show dense gas in their central regions (e.g., Henkel, Baan, & 
Mauersbcrgcr 1991), so it is natural to ask if one could apply the above analysis and obtain 
the resulting steady-state density of the dense, virialized clouds. CaroUo and Stiavelli (1998) 
have observed using HST the nuclear surface brightness profile of a large sample of spiral 
galaxies. A considerable fraction of these spirals have a fiat core luminosity profile close to 
their center. The galactic tidal field for these galaxies can become compressive within the 
core. We plan to study these galaxies in further detail in a forthcoming paper. In particular 
we would like to examine how the molecular cloud densities in these galaxies compare 
with those predicted by our work. This may help us decide whether the so-called diffuse, 
inter-cloud medium in normal and starburst galaxies (e.g., Bally et al. 1988; Jog k. Das 
1992) actually consists of bound, virialized clouds. 

2. However, the mass distribution and hence the dynamics in the centres of spirals is more 
complex due to the common presence of bars and triaxial bulges. We find that for the 
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Milky Way Galaxy, assuming the bar potential as in Binney et al. (1991), the field is not 
compressive. A typical bulge potential (e.g., Sellwood & Wilkinson 1993) also does not result 
in compressive fields. Physically, this could be because a bar is a strong deviation from an 
axisymmetric potential, and the field near the bar is not smooth. In fact, for a molecular 
cloud close to a strong bar, the galactic tidal field is shown to cause an internal heating of 
the cloud (Das & Jog 1995). We will discuss the barred case in a future paper. 

3. We have only considered the radial tidal field in this paper. The tidal fields along the other 
two orthogonal directions are always compressive (see Section 2.1), and have magnitudes 
comparable to the radial case. We are justified in neglecting this asymmetry since the gas 
clouds are compressible and would tend to have an isotropic velocity dispersion, and hence 
the net eff'ect of the compressive fields would be seen in an average sense. In the case of the 
ultraluminous galaxies, this asymmetry would be particularly stronger since the height of 
the disk is small and hence using a 3-D distribution is not strictly correct. 

4. In the centers of ultraluminous galaxies, the gas constitutes a significant fraction of the mass 
(> 25%) and hence the gravitational potential. Hence a gas cloud does not respond to the 
stellar potential alone. Since we have used the observed gas rotation curve in obtaining the 
potential, we have used the potential as seen by a gas cloud. 

6. CONCLUSIONS 

We have investigated a simple but novel physical process which can compress gravitationally 
bound, virialized molecular clouds in the centers of galaxies. The main conclusions of this study 
are the following. 

1. The galactic tidal field changes from being disruptive to compressive in the centers of fiat 
core early type galaxies and ultraluminous galaxies. 

2. The compressive tidal force and the self gravity of the cloud are balanced by the internal 
pressure of the cloud. We have used the virial theorem to derive a simple expression for the 
minimum mean density of a gravitationally bound cloud in the compressive core of a galaxy. 
The additional inward/compressive field increases the velocity dispersion required to satisfy 
the virial equilibrium, and hence the critical density required for stability in presence of a 
compressive force is higher for a gravitationally bound, virialized cloud. For lower densities, 
the gas within that region will be diffuse and gravitationally unbound. Such gas will not 
lead to star formation. 

3. We have applied our results to a sample of flat core, early type galaxies taken from Faber 
et al. (1997). The gas densities predicted by the model varies with core size and cloud 
position in the center of the galaxy, and the typical density range observed within the central 
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lOOpc is (n) ^ 10^ — lO^cm"^. The galaxies in our sample were chosen so that they have 
an observational evidence of molecular gas in their centers. However, the spatial resolution 
of these data are not yet adequate and only give a beam-averaged lower limit on the gas 
densities. Future high resolution studies of such galaxies are highly recommended. 

4 We have also applied our results to the molecular gas disks observed in the centers of 
ultraluminous galaxies (Downcs & Solomon, 1998). From their observed rotation curves, we 
show that they have a compressive core within their central region. The minimum mean 
molecular cloud densities predicted from our model are in the range of lO^cm"^ to lO^cm"^ 
and this range matches well with the observed gas densities in these galaxies. Also, if the 
tidal field is compressive and not disruptive as previously thought in the literature, this 
molecular gas can exist in the form of bound clouds. Hence, these clouds arc closer to star 
formation than unbound, diffuse molecular gas. This strengthens the scenario proposed by 
Downes &; Solomon (1998) for starbursts powering the ultraluminous galaxies. 

M.Das would like to thank R.Nityananda for useful discussions during the course of this 
work, and we would like to thank M. Valluri and S. Faber for very useful e-mail correspondence 
regarding this work. We thank the anonymous referee for helpful comments which improved the 
comparison with observations of ultraluminous galaxies. We thank S. Sridhar for comments on 
the choice of reference frame which we have addressed in the Appendix. 
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A. APPENDIX A. CLOUD EQUILIBRIUM IN THE INERTIAL AND 

ROTATING FRAMES 

In this appendix we derive the condition of cloud stabihty (see eq.[10]) in the galactic reference 
frame and in the rotating frame moving with the cloud. We then discuss which equation was 
used in the paper and the reasons behind choosing that particular form. The coordinate of the 
cloud with respect to the center of the galaxy is Rq. The i'*'^ clump in the cloud has a position 
coordinate Ri wrt the galaxy center and ri wrt the cloud center, so that Ri = Rq + rj. The 
number of clumps in the cloud is n so that the index i = 1 to n. 

(i) Galactic Frame : The origin of the coordinate system is the center of the galaxy. For the i'**^ 
clump within the cloud the total velocity Vjg is given by, 

Vig = Vo + Vi (Al) 

where Vq is the rotational velocity of the cloud about the galaxy center and Vi is the random 
motion of the i^**^ clump. The kinetic energy of the cloud K, is given by, 

n / " \ 

2K = MM^ + + 2Vo. 5] Pi (A2) 

i=l \i=l ) 

But when the the last term is averaged over time, the random velocities of the clumps makes the 
term average out to zero. Hence equation(A2) reduces to, 

n 

%K) = McFo' + (E "^i^i") (A3) 

i=l 

To determine the right hand side of equation (8), we first determine the force on the i^**^ clump, 
which is the sum of the galactic potential and the self gravity of the cloud. 

(A4, 

where rjj = (ri — rj). Expanding the force term about the radius Rq gives the tidal term so 

that EiLi^i.ri reduces to, 

E Fi.ri = -^Ri. E m,ri + To ^ rn,r^ - ^^=' (A5) 

Using the relation McRq = Y^i=i ^i^h the equation of virial equilibrium for the cloud reduces to, 

{j^m.v,'') = ^^=15^=1 - n{±m,r^) (A6) 

We have used this form of the virial equation (see eq.[10]) to determine the minimum cloud density 
in a compressive tidal field. 
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(ii) Rotating frame : In this case the coordinate system is centered in the cloud so that the i'*^ 
clump has a position vector ri. The force in the rotating frame is given by, 

Fi = -mi [ Ri - 2L/ Tmr-^ij " 2mi(a;o x vh) - rriiioo x {loq x Ri) (A7) 



j=l I "^ij 



where loq is rotation velocity of the cloud about the galaxy cente r i.e. I "^o |= ^ and Vri is the 
random motion of the i'*^ clump in the rotating frame. The third term on the right hand side of the 
equation is the Coriolis force and the fourth term is the centrifugal force. We expand the force term 
(^1^^ about Ro and put Ri = Ro + ri. Since (|h)o ~ '*^o^Ro and ojq x (wq x Rq) = — wq^RoRo; 
the force Fi becomes, 

Fi = TomiTi - ^ - — ^-/rij - 2mi{uio x vh) - niiUJo x (cjq x (A8) 



j=l I '^ij 



Hence, Fi.ri is given by. 



Fi.ri = Torrurf ^ - — ^rij.ri - 2mi{uo x Vri).ri - ruiUo x (a;o x ri).ri (A9) 

i=i I ^ii I 

Now, ujQ X {ooq X ri.)ri = — (wq x rj)^ and {ooq x Vri).ri = — a;o.(ri x Vri). So equation (A9) 
becomes, 

Fi.ri = Toniirf - ^ - — ^rij.ri + 2uo.{miri x v^) + mi{LOo x rif (AlO) 
j=i I fij I 

Also, Li = mi(ri x Vi) is the angular momentum of the clump measured in the inertial or galactic 
frame and Lir = mi(ri x Vri) is the cloud angular momentum as measured in the rotating frame. 
The two are related through the following relations (e.g. Goldstein, 1980), 

Vig = Vri + Wo X Ri (All) 

Since Vig = Vq + Vi, Ri = Rq + ri and Vq = loq x Rq, we obtain the relation, 

wo-Li = wo-Lir + rrii^LOo x n)^ (A12) 

However, observations indicate that clouds are practically non-rotating, so the net angular 
momentum of a cloud as observed in the inertial frame is ~ (Blitz, 1991). Applying this 

to equation (AlO) we get the following relation, 

E(^"') = ^^'P=' - To f: m.n' + ^ m,{uo x ri)^ (A13) 

So the equation of equilibrium is different from that in the inertial or galactic frame by one term 
(see eq.[A6]). However, using the relation vh^ = (vi — loq x ri)^, we obtain, 

E rriiVri^ = E miVi"^ + ^ mi{uio x n) (A14) 
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which when substituted in equation (A13) leads to an identical relation for cloud equilibrium as 
obtained in the galactic frame and as expected physically. 

We have chosen to use the Galactic or inertial frame equation eq.[A6] and not eq.[A13] simply 

because the velocity is observed from the galactic reference frame and not from within the cloud 
(which is the rotating frame). Hence we use the velocity dispersion (X^jVi^) and not (Z^iVn^). 
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TABLE 1 



Model Cloud Densities for Early Type Galaxies 



Galaxy 


Galaxy 


Core/Power 


Break Radius 


Galaxy Mass 


Predicted H2 Number 


Name 


Type 


Law Type 


Ro (pc) 


(Mo) 


Density (cm^) at R'=0.4 


NGC 4594 


Sa/Sb 


Power Law 


4 


1.1 X 10" 




NGC 4697 


E6 


Power Law 


240 


1.5 X 10" 




NGC 4472 


E/SO 


Flat Core 


178 


8.4 X 10" 


lO"* 


NGC 4649 


El-2 


Flat Core 


263 


9.9 X 10" 


5.10^ 


NGC 1400 


El/so 


Flat Core 


35 


2.4 X 10" 


5.10« 


NGC 1316 


SO 


Flat Core 


36 


2.9 X 10" 


5.10^ 


NGC 1600 


E 


Flat Core 


759 


1.5 X 10^2 


300 


NGC 3379 


E 


Flat Core 


83 


9.8 X 10^° 


10"* 


NGC 4486 


E 


Flat Core 


562 


9 X 10" 


420 



1 
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TABLE 2 



Model Cloud Densities for Ultraluminous Galaxies 



Galaxy 
Name 


Rotation Curve 
Turnover Radius (pc) 


Rotation Curve 
Turnover Velocity (kms~^) 


Predicted H2 Number 
Density (cm^) 


00057+4021 


240 


250 


7 X 10^ 


02483+4302 


270 


270 


6 X 10^ 


VIIZw31 


290 


290 


6 X 10^ 


10565+2448 


230 


220 


6 X 10^ 


Mrk231 


75 


345 


10^ 


Arpl93-disk 


220 


230 


7 X 10^ 


Mkr273-disk 


70 


280 


10^ 


Arp220-disk 


200 


330 


2 X 10^ 


Arp220-west 


35 


300 


5 X 10^ 


Arp220-east 


50 


350 


3 X 10^ 


17208-0014 


310 


260 


5 X 10^ 


23365+3604 


340 


260 


4 X 10^ 



1 



